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ABSTRACT 

Recent observations of blue stragglers by iDe Marco et al. (200 li have revealed continuum deficits 
on the blue side of the Balmer discontinuity, leading these authors to infer the presence of discs around 
the stars. This intriguing possibility may throw light on aspects of the mechanisms responsible for 
at least some of these objects; current theories of blue straggler formation invoke stellar collisions or 
interacting binaries, both of which appear capable of forming a circumstellar disc. 

However, by synthesizing photospheric spectra for models of rotating blue stragglers, we demonstrate 
that the Balmer jump enhancements can be wholly attributed to the influence of oblateness and gravity 
darkening on the formation of the continuum. Therefore, we are led to conclude that the observations 
of De Marco et al. can be ascribed a more prosaic explanation, that of rapid stellar rotation arising 
from the merger/interaction formation process. 

Subject headings: blue stragglers - stars: rotation - stars: atmospheres - techniques: spectroscopic 


1. INTRODUCTION 

Blue stragglers (BSs) are cluster stars having anoma¬ 
lous evolutionary histories. T hey are blue stars of in¬ 
termediate mass (a few M™. lSha,ra, et all I1997 L with 
typical surface te mperatures of ~ 6, 000-10, 000 K (e.g., 
iDeng et "all 11999th Application of single-star evolution 
theory indicates an age too young to be expl ained b y the 
age of their parent cluster (see the review bv lLiviolll99,4 
and references therein). This apparent paradox has been 
explained by two production paths for BSs: (i) collisions 
be tween two l ower -mass stars (e.g. see the simulations 
bv lSills et al.l200llb and (ii) mass transfer in moderately 
wide binary stars (e.g.. iBellazzini et alJ l20021. Both of 
these processes are likely to cause the remnant BS to spin 
up. Population studies of BSs have shown that both of 
these rou tes may be necess ary to account for the obser¬ 
vations llDavies e t al.ll2004D 

Recently. iPe Ma~m"et~al3 ( 2004 . hereinafter DM04) 
have presented Hubble Space Telescope ( HST ) observa¬ 
tions of BS spectra, and have argued that they detect the 
signature of a circumstellar disc around some of their tar¬ 
gets. Disc formation is a common theme in BS generation 
and hence the DM04 observations could shed light on the 
different evolutionary processes responsible for these ob¬ 
jects. The evidence advanced for the presence of discs 
rests largely on an apparent continuum deficit (by ~ 5- 
10%), falling on the short wavelength side of the Balmer 
discontinuity; DM04 interpret this deficit - and the cor¬ 
responding enhancement in the magnitude of the Balmer 
jump - as arising from circumstellar absorption by the 
presumed disc. 

On account of their putative formation mechanisms, 
BS are expected to exhibit moderate to high rotation 
rates. In this letter we investigate the effect of such rapid 
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rotation on the photospheric continuum, by developing 
a spectral synthesis model that correctly accounts for 
the oblateness and gravity darkening arising from the 
centrifugal force. We describe the model in and use 
it in to explore the evolution of the Balmer jump as 
the rotation rate is varied. In 0we then apply the model 
to the blue straggler M3-17, demonstrating that it can 
successfully reproduce the observations by DM04. We 
discuss and summarize our findings in m 


2. THE SPECTRAL SYNTHESIS MODEL 

In a non-rotating, spherical star, the atmosphere across 
the entire surface may be characterized by a single value 
each for the effective temperature T e g and gravity g. 
With the introduction of rotation, however, the out¬ 
ward pull of the centrifugal force distorts the star into 
an oblate spheroid. Across the surface of this spheroid, 
the effective gravity g e g (composed of the vector sum of 
the Newtonian gravity and the centrifugal acceleration) 
is non-uniform, decreasing toward the equatorial regions 
where the centrifugal force is strongest. Like wise , i n ac- 
cordance with the gravity darkening law of Ivon Zeinell 
lTf92l . the effective temperature decreases toward the 

l 

equator in accordance with the relation T e g oc <? e 4 ff . 

To model such an oblate, gravity-darkened star, we set 
up a raster grid that divides the surface 1 into c. 14,000 
pixels 2 , each covering an area 0.015 R p square, where R p 
is the star’s polar radius. We calculate the local effective 
gravity associated with the surface point to which each 
pixel corresponds, and we likewise assign a local effective 

1 Defined by an isosurface of_tdie_effectiver>otentiaI within the 
Roche approximation (see, e.p.. ICranmertlll 996K 

2 The exact number of pixels composing the surface varies with 
the degree of centrifugal distortion of the star. 
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Fig. 1. — Synthetic spectra from rotating stars: the top panel 
shows a non-rotating model, and (decreasing downward for wave¬ 
lengths less than 3700 A) models for n e q/n cr it = 0.3 and 0.5, while 
the lower panel have models for Heq/r’crit = 0.7 and 0.9 (increasing 
upward). All models are equator-on (inclination i = 90°). The 
insets are the models presented over a smaller wavelength range to 
highlight the flux shortward of the Balmer Jump. The dotted line 
on all panels shows the spectrum of M3-17. 
temperature 
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foil owing Ivon Zeinef s law. The overall normalization of 
these temperature data is specified by the notional ef¬ 
fective temperature T^o that the star would have if 
it were non-rotating, under the (reasonably-accurate) 
ansatz that the stellar bolometric luminosity remains in¬ 
variant as the rotation rate changes. The symbol £i de¬ 
notes the surface-area weighted gravity of the distorted 
star, which depends amongst ot her thing s on the rotation 
rate (see equations 4.22-4.24 of[C ranme il l 19961) . 

Knowing the effective temperature and gravity of each 
pixel, and furthermore the local projection cosine of the 
surface normal onto the line of sight, we interpolate the 
observer-directed emergent flux in a T e ff-log g-g grid of 
precomputed angle-dependent intensity spectra. By co¬ 
adding the flux data from all pixels, weighted by their 
projected area (0.015i? p ) 2 and Doppler shifted by the 
line-of-sight velocity due to rotation, we thereby build 
up a disc-integrated spectrum for the entire distorted, 
gravity-darkened star. For the calculations presented in 
the following section, we adopt an intensity spectrum 
grid calculated using the SYNSPEC spectral synthesis code 
of I. Hubeny and T. Lanz. The spectra incorporate 
lines due to H, He, C, N, O, Si, Mg, and Ne, and are 
based on the am20ak2-odfnew grid of line-b lanketed LTE 
model atmospheres published by iKuruc 2 (these 

have H/He = 0.34 by mass and are alpha enhanced by 
0.4dex and metal depleted by 2.0dex). 


3. INFLUENCE OF ROTATION ON THE BALMER JUMP 
AND COLORS 

In Fig. n we present synthetic spectra covering the 
wavelength range 3,200-4,500 A, calculated for equa¬ 
torial rotation rates u e q/u cr it = 0.0,0.3, 0.5, 0.7 and 


Fig. 2.— The Balmer jump magnitude Ab plotted as a function 
of effective temperature at fixed log# = 4.0 (left), and as a function 
of gravity at fixed T e ff = 10, 000 K (right). 

0.9 using the model we describe above; here, u cr it = 
v / 2GM*/3i? p is the critical rotation velocity at which 
the equatorial centrifugal force balances gravity. The 
underlying star has a mass M* = 1.35 Mq, polar radius 
R p = 2AR & and ‘non-rotating’ effective temperature 
T e ff,o = 10 , 000 K, these parameters being chosen to co¬ 
incide with those given by DM04 for the star M3-17. 
The model star is viewed equator-on, and we normalize 
all spectra to have a unit flux at 4,200 A. 

It is clear that the synthetic spectra are all similar 
longward of the Balmer discontinuity at 3,647 A (see be¬ 
low) but that the magnitude of the jump - which we 
characterize throughout via the ratio Ab = / 4 , 20 o// 3,630 
between the continuum fluxes at 4,200 A and 3,630 A 3 - 
varies noticeably with the rotation rate. In interpreting 
this behaviour, we recall that the spectrum of a rotating 
star is a composite, made up from contributions cover¬ 
ing a range of effective temperatures and gravities. To 
illustrate the sensitivity of the Balmer jump against such 
variation in temperature and gravity, Fig. |2] plots Ab as 
a function of both T e g and log < 7 , for a spectrum syn¬ 
thesized from a non-rotating plane-parallel atmosphere 
model. 

The temperature dependence of Ab (left-hand panel) 
exhibits a sharp peak around T e s ~ 8 , 500 K; at temper¬ 
atures cooler than this turnover, the appearance of H~ 
bound-free opacity - which preferentially absorbs con¬ 
tinuum photons toward longer wavelengths in the opti¬ 
cal and UV (see, e.g. Fig. 8.3 of lGravlll992h - tends to 
suppress the flux redward of the Balmer discontinuity, re¬ 
sulting in a reduction of Ab- Likewise, at temperatures 
hotter than 8,500 K, the progressive depopulation of the 
n = 2 (Balmer ground-state) level of neutral hydrogen 
removes bound-free continuum opacity blueward of the 
discontinuity, again resulting in a reduction of Ab- 

Similar behaviour is exhibited in the gravity depen¬ 
dence of the Balmer jump (right-hand panel of Fig. 0 ; 
H~ number densities are enhanced at high gravities, due 
to the corresponding increase in electron density, while 

3 These points being chosen as well separated from spectral lines 
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Fig. 3.— The Balmer jump magnitude Ab plotted as a function 
of rotation rate Ueq/u C rit 5 f° r the stellar model representative of 
M3-17 (see text). Data for different inclinations are shown; the 
horizontal dashed line indicates the jump magnitude observed in 
M3-17. 

n = 2 level populations are depleted at low gravities, due 
to the reduction in collisional recombinations which re¬ 
plenish these populations. Together, these processes are 
responsible for the decline in Ab toward both low and 
high g. 

In Fig. EH we illustrate how the combined T e g/\ogg sen¬ 
sitivities of the Balmer-jump magnitude come together 
in a rotating star, by plotting Ab as a function of rota¬ 
tion rate. We use the same stellar parameters as previ¬ 
ously, but in the present case show data for both pole- 
and equator-on aspects. In each, Ab increases steadily 
up to a rate v eq /v cr it ~ 0.77, due to the reduction in 
equatorial effective temperature and gravity from their 
non-rotating values T e s = 10, 000 K, logoff = 3.81. At 
this point, the equatorial reduction of T e g and g e s is so 
pronounced that the low temperature/gravity regime of 
Ab is reached, and the magnitude of the Balmer jump 
then decreases rapidly toward even-higher v eq . However, 
in every case, Ab for the rotating models is larger than 
in the non-rotating limit. 

How then does rotation change the flux distribution 
redward of the Balmer Jump? In Fig. we present a 
color-color diagram of the rotating models. The fluxes in 
the three photometric bands (centered at 3660A, 4200A 
and 5450A) are averaged over 120A, 80A and 120A 
widths resp ective ly, and we use the zero -point magni¬ 
tudes from iBessell Castelli fc Pled (11998H - this is ex¬ 
actly the same as the procedure used DM04 to allow a 
direct comparison to be made. A grid of non-rotating 
model atmospheres is plotted (dotted lines) for temper¬ 
atures between 8,500K-12,000K, and for logg=2. 5-5.0. 
Also we have calculated pole-on and equator-on rotat¬ 
ing atmosphere models (for v eq /v CI it = 0 . 0 - 0 .9 in steps 
of 0.1) for the [4200]-[5450] extreme blue edge of the 
non-rotating grid. The equator-on models (filled cir¬ 
cles) sweep down and to the right (almost appearing 
like models becoming cooler at constant logg), with the 
lower-temperature models eventually sweeping upward 
in a broad “u” shape. The loci of the pole-on models 



Fig. 4.— Color-color figure for both rotating and non-rotating 
atmosphere models. The dotted lines (with effective temperatures 
and gravities marked) correspond to models with zero rotation. 
The solid lines correspond to models with rotation rates of v/v CT [t = 
0.0-0.9 in steps of 0.1 for pole-on (squares) and equator-on (circles) 
models. 

(filled squares) initially move vertically down with in- 
creasing rotation (i.e. no change in the [4200]-[5450]) 
with the cooler models then reversing their track in a 
“v” shape. The maximum change in [3660]-[4200] be¬ 
comes more pronounced for higher temperature models, 
reaching 0.06mag for the 12,000K model. 

Two conclusions can be drawn from this figure: first, 
that for pole-on models, rotation can cause the Balmer 
jump (signified by the [3660]-[4200] color) to be too large 
whilst keeping the [4200]-[5450] color constant; second, 
if they are rotating, stars can exist with colors which are 
apparently inconsistent with non-rotating model atmo¬ 
spheres (i.e. the lower left-hand extremum of the non¬ 
rotating atmosphere grid is extended a little). Finally we 
note that the position of a star in the color-color plane 
does not yield a unique pair of effective temperature and 
gravity values - two color measurements cannot produce 
three unique values of effective temperature, gravity, and 
rotation. 

4. APPLICATION TO M3-17 

The observed projected rotation velocity of M3-17, 
usini = 200 ± 50kms _1 , along with the stellar param¬ 
eters furnished by DM04, implies that this object is a 
rapid rotator, having v eq sin i/v CT it = 0.75 ± 0.18. The 
Balmer jump magnitude derived from the HST spectrum 
of this object is Ab = 2.47, a value that we indicate in 
Fig. 0 by the horizontal dashed line. Clearly, there are 
multiple rotating models that can fit the observations, 
without the need to invoke a circumstellar disc; for ex¬ 
ample, a match to Ab can be achieved with an edge-on 
(i = 90°) model rotating at Weq/'Ccrit ~ 0.4 or 0.9, and 
other solutions can be found at lesser inclinations. 

We further illustrate this point in Fig.0 where we plot 
the archival spectrum of M3-17 over our model data. The 
non-rotating synthetic spectrum clearly has a surfeit of 
continuum flux shortward of the Balmer discontinuity, 
which led DM04 to diagnose the presence of a circum¬ 
stellar disc. However, our equator-on rotating models 
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at Ueq/^crit = 0.4 and 0.9 clearly provide a good fit to 
the continuum flux, suggesting that rotation alone may 
suffice to explain the Balmer jump. In support of this 
conclusion, we note that the rotational reddening of the 
Paschen continuum, as seen for the equator-on curves in 
Fig. El actually enables the rotating models to fit the ob¬ 
served [4200]-[5450] color of M3-17 better than the non¬ 
rotating model of the same temperature. We do observe 
that our models are unable to reproduce the strengths 
of the Balmer-series absorption lines seen redward of 
3,647A; this problem reflects the fact that we specify the 
model temperature a priori from that of DM04. The is¬ 
sue can be resolved by allowing T^o (and other parame¬ 
ters) to vary under the control of a suitably-chosen fitting 
statistic. However, we do not undertake such fine tun¬ 
ing here, because we are not attempting to derive exact 
atmospheric parameters for M3-17 here. Rather, our ob¬ 
jective is to highlight that high rotation has a significant 
effect on the continuum, and may mimic the presence of 
an obscuring disk. 

5. DISCUSSION & SUMMARY 

Using a spectral synthesis approach that correctly 
treats oblateness and gravity darkening, we have demon¬ 
strated that apparent flux deficits shortward of the 
Balmer discontinuity can be attributed wholly to rapid 
rotation. This casts doub t on t he disc hypothesis ad¬ 
vanced bv lDe Marco et ah ( 12004) to explain the Balmer 
jump anomalies seen in the spectra of BSs. At least in the 
case of rapid rotator M3-17 (usini = 200 ± 50kms _1 ), 
the continuum is well fit by our model, leaving little rea¬ 
son to invoke the presence of a disc. 

4 By this, we mean models that do not include oblateness or 
gravity darkening. 


In addition to M3-17, five other stars (out of a total 
of 50) in the DM04 sample show Balmer jumps that are 
enhanced with respect to non-rotating models 4 . Of these 
five, only NGC6751-11 has a measured (upper limit) pro¬ 
jected rotation velocity, v sin i < 50kms _1 . Due to the 
unknown projection factor sinz, the upper limit on the 
intrinsic equatorial velocity 'C e q may be much larger, and 
it is entirely possible that, along with the remaining ob¬ 
jects with anomalous Balmer jumps, this star is a rapid 
rotator. 

Accordingly, Balmer jump anomalies seen in BSs may 
all be attributable to rapid rotation. This hypothesis 
is lent some support by the fact that the proposed for¬ 
mation mechanisms for BSs - collisions or mass transfer 

both involve the deposition of significant amounts of 
angular momentum on the star , thereby spinning it up 
(see, e.g.. lPols fc Marinu^ll994l) . 

On a final note, we stress that we do not (and can¬ 
not) rule out the possibility of discs around BSs. In¬ 
deed, observations of emission lines in these objects, 
bear ing a close resem blance to those seen in Be stars 
fsee iMermilliodj ll982 ). indicate that there almost cer¬ 
tainly are circumstellar discs around some of them. How¬ 
ever, the present paper demonstrated why Balmer jump 
anomalies cannot be used as a clear and unambiguous 
diagnostic for BS discs. 
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